Capacitively Enhanced Thermal Escape in Underdamped Josephson Junctions 
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We have studied experimentally the escape dynamics in underdamped capacitively shunted and 
unshunted Josephson junctions with submicroampere critical currents below 0.5 K temperatures. 
In the shunted junctions, thermal activation process was preserved up to the highest temperature 
where the escape in the unshunted junctions exhibits the phase diffusion. Our observations in the 
shunted junctions are in good agreement with the standard thermal activation escape, unlike the 
results in the unshunted junctions. 



Underdamped Josephson junctions are routinely used 
as threshold current detectors to readout quantum infor- 
mation in a superconducting circuit by measuring switch- 
ing events from the superconducting to the normal state 
(see, e.g. [1]). In addition, this technique has been used 
to identify between geometric and dynamic phases in a 
superconducting charge pump 0, and also to measure 
collective quantum phase slips in the ground state of a 
Josephson junction chain The detector sensitivity is 
important in the measurement. One can improve it by 
reducing the critical current of the junction, but then the 
system typically enters the phase diffusion regime, which 
is not desirable since in this context one usually wants to 
fix the phase difference of the superconducting circuit. In 
this letter, we present a technique to prevent the phase 
diffusion by adding a large shunt capacitance with the 
help of atomic-layer deposition (ALD). 

The Josephson junction dynamics is well described by 
a model of a phase particle in a tilted washboard poten- 
tial [1]. In this model, as the bias current / is increased 
toward the critical value Ic^ the barrier height AU is 
lowered and the particle, oscillating with the plasma (an- 
gular) frequency ujp = {2elc /Cjh)^^'^ in a potential well, 
where Cj is the junction capacitance, can escape out of 
the well by either thermal activation (TA) or macroscopic 
quantum tunneling (MQT) process. The crossover be- 
tween MQT and TA regimes occurs at the temperature 
Above Tcr, the dominant escape 
mechanism is TA. The junction behavior also depends 
on the magnitude of the quality factor Q of the junction 
at its plasma frequency. In the case of junctions with 
Q <C 1, called overdamped, the escape dynamics occurs 
over a dissipation barrier [6|, which is undesirable for a 
threshold current detector. For Q 1, so-called under- 
damped case, the phase particle moves down freely or dif- 
fusively from one well to another depending on the effect 
of dissipation. When dissipation is weak, the free running 
particle is still governed by TA process, while when dis- 
sipation is strong, the escape dynamics enters the phase 
diffusion regime [Bl, 0-Q • The appearance of phase dif- 
fusion in a hysteretic junction was studied extensively 
two decades ago by Martinis and Kautz [ll|, [l2[ . How 
can one avoid the underdamped phase diffusion regime 
and maintain TA process at high temperatures? Accord- 
ing to the phase diagram in Ref. [sj, decreasing Tcr can 



help the junction escape from the phase diffusion regime 
by enhancing TA. This idea can be realized by either 
decreasing Ic or increasing Cj. Among these two solu- 
tions, however, only the second one is applicable, since 
the junction with small Ic has also a low Josephson cou- 
pling energy ^j, which lowers /SXJ . Low /SXJ implies a 
small tilt of the potential at the point of escape, hence 
the particle gains a small amount of energy when es- 
caping and can get trapped in the next well instead of 
running down the potential. Accordingly, Cj is the only 
parameter remaining for decreasing Tcr- For a typical 
junction, the capacitance per unit area is 45 fF//im^. 
With this value, although one can implement a large Cj 
by increasing the junction area, it becomes impractical 
for the detector junction, since Ic is also proportional to 
the area. Alternatively, a large shunt capacitance can be 
added in parallel to Cj. 

Motivated by this idea, we implemented a shunt capac- 
itance Csh by growing AI2O3 layer of 5.5 nm thickness d 
using ALD on top of an Al ground plane of 40 nm thick- 
ness. Using electron beam lithography and shadow angle 
evaporation [lo|, Al/AlO^c/Al Josephson junctions were 
fabricated on the ALD oxide. The insulating AlO^c layers 
for the tunnel junctions were formed by standard room 
temperature oxidation between the evaporation steps of 
the two Al layers, which allows us to manipulate the junc- 
tion resistance and Ic by varying the oxidation pressure. 
For large capacitance, two rectangles with 30x15 /im^ 
area were defined in the leads near the junctions, over- 
lapping with the ground plane, as shown in Fig. 1(b). 
The total overlap area Atot between the two parts con- 
nected to the junction on the top and the ground plane 
at the bottom is roughly 500 jim^ each. As a result, the 
Josephson junction is capacitively shunted by the paral- 
lel plate capacitor which consists of two shunt capacitors 
Co in series and on each side of the junction. The total 
shunt capacitance is C^h = Cq/2 where Co = eAioi/d 
and e is the permittivity of AI2O3. A similar tech- 
nique for adding shunt capacitance to proximity Joseph- 
son junctions was reported in Refs.0,1^, where the ca- 
pacitor changes the damping strength in superconductor- 
two-dimensional electron gas-superconductor Josephson 
junctions. The leakage resistance from the junction to 
ground plane was larger than 10 G^l. Both the unshunted 
and shunted junctions were fabricated at the same time 
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FIG. 1: Scanning electron micrographs of (a) an unshunted 
and (b) a capacitively shunted junction. In (b), the junc- 
tion and leads (purple, in foreground) overlap with Al ground 
plane (encircled by solid white line, in background). The sizes 
of the junctions in the panel (a) and (b) are 0.4 x 2.5 /im^ 
and 0.4 x 1.5 /xm^, and the circuit schematics of the junctions 
correspond to (c) and (d), respectively. 



FIG. 2: Measured switching histograms of (a) the unshunted 
junction Al and (b) the shunted junction A2. The bath tem- 
perature Tbath was set to 58, 120, 160, 216, 270, 321, 377, and 
430 mK for the different sets of data from the right to the 
left. 



on the same chip. Note that the total capacitance Ctot is 
Ctot — Cj in the unshunted junction and Ctot — Cj^C^h 
in the shunted junction, as shown in Fig.l (c) and Fig. 
1 (d), respectively. All measured data were taken in a 
^He-^He dilution refrigerator with a base temperature of 
50 mK. 



TABLE I: Parameters of the samples. Rn is the normal- state 
resistance, Ic is the critical current, and Ej is the Joseph- 
son energy obtained from Ambegaokar-Baratoff formula [ij] . 
Ctot is the total capacitance including the junction capaci- 
tance Cj and the shunted capacitance Csh- 



Sample 


Rn Ic 

Device . . . 

(Q) (nA) 


Ej 

(K) 


Ctot 

(pF) 


A 


Unshunted (Al) 612 
Shunted (A2) 828 


513 
379 


12.2 0.045 
9.0 3.2 


B 


Unshunted (Bl) 775 
Shunted (B2) 746 


405 
421 


9.7 
10.0 


0.045 
3.2 



In Table I, we present the characteristics of the de- 
vices studied here. The normal-state resistance Rn was 
obtained from current-voltage (IV) measurement, which 
yields In and Ej according to Ambegaokar-Baratoff for- 
mula [13|- By the standard parallel plate capacitance 



calculation based upon the information above and the 
relative dielectric constant 7.8 of AI2O3 in ALD process- 
ing, Csh was estimated to be 3.14 pF. We also estimated 
Cj from the specific capacitance of 45 £F/(/im)^ and the 
junction areas given in Figs. 1(a) and (b). Consequently, 
the total capacitances of the shunted and unshunted junc- 
tion were 3.2 pF and 45 £F, respectively. 

We investigate the escape dynamics by applying a set 
of current pulses with a fixed amplitude and duration 
through the junctions and by determining the probability 
for the devices to switch from the zero- voltage state to the 
finite- volt age state. At each amplitude 200 current pulses 
were applied. The escape probability P was measured as 
a function of the current pulse amplitude /, which yields 
cumulative histograms of switching current. Figure 2 
shows the measured switching histograms of both the un- 
shunted (Al) and the shunted junctions (A2) in Sample 
A at different bath temperatures. A similar measurement 
was carried out with Sample B, but the histogram curves 
are not shown here. The histogram is typically charac- 
terized by the switching position 750% = I{P = 0.5), and 
the width Al = I{P = 0.9) - I{P = 0.1), the measure- 
ment results of which are shown in Fig. 3. Note that at 
the lowest temperature, the switching position of A2 is 
lower than that of Al, as shown in Fig. 3(a), and Sample 
B behaves in an opposite manner as shown in Fig. 3(b). 
We observed that on increasing bath temperature Tbath 
the mean position shifts down in all cases, whereas the 
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FIG. 3: The temperature dependent variation of the mean po- 
sition (/5o%) in (a)-(b) and the width (A/) of the histograms 
in (c)-(d). The experimental data in panels (a) and (c) and in 
(b) and (d) were obtained from Sample A and Sample B, re- 
spectively. In each graph, diamond (o), circle (o), and dashed 
line correspond to the experimental data of the unshunted 
junction (Al and Bl), the experimental data of the shunted 
junction (A2 and B2), and the theoretical simulation of the 
shunted junction by the TA model, respectively. The simu- 
lation was performed with Ic of 379 nA for Sample A and 
with 369 nA for Sample B as fit parameters, and Ctot = 3.17 
pF was employed. The latter Ic differs by 12 % from the 
Ambegaokar-Baratoff value. In the TA model, AI does not 
depend on Ctot, or at least the dependence should be very 
week. 

temperature dependence of the width is very different de- 
pending on whether the junction is capacitively shunted 
to the ground plane or not. According to the standard 
MQT and TA discussed above, Tcr - 11 mK <C Tbath 
in the shunted junctions (A2 and B2). Thus the role of 
the MQT is negligible and instead, TA dominates the es- 
cape dynamics. In the TA model, the width obeys AI 

2 /3 

oc T^ath- Consistent with the prediction of the model, 
we observed that in the shunted junctions the widths 
are increasing as a function of temperature. The dashed 
lines in Fig. 3 show that the theory is in quantitative 
agreement with our observations. For the corresponding 
unshunted junctions (Al and Bl), on the other hand, on 
increasing Tbath the widths are more or less constant up 



to Tcr ^ 200 mK, as predicted with the MQT model, but 
they decrease at high temperatures. It was reported that 
the decreasing AI results from the influence of dissipa- 
tion by a retrapping process [5|,0-[2|- This indicates phase 
diffusion. The TA model cannot explain the phenomenon 
because it neglects the effect of dissipation. Several mod- 
els have been proposed to account for the phase diffusion 
of underdamped Josephson junctions incorporating the 
retrapping process in the TA model [El. ItMoI.!!^ . 

Kivioja et al. [5] claim that when / is below a maxi- 
mum possible phase diffusion current Im = 4:1c /ttQ^ the 
phase particle is retrapped after escape. The observed 
Im is 265 nA in Al and 190 nA in Bl, obtained from 
/50% curves at a certain temperature where AI starts to 
decrease in Fig. 3. With those values and Ic given in 
Table I, Q of the unshunted junctions is estimated to be 
2.5 in Al in and 2.7 in Bl. Note that this Q may differ 
from that introduced earlier since here it corresponds to 
low frequency dynamics of the junction. For the shunted 
junctions, however, it is not simple to determine Q with 
the present information. Nevertheless, we believe that Q 
must be at least an order of magnitude larger than those 
of the unshunted junctions, because of the large Ctot- 

In conclusion, we developed a technique to add a 
large shunt capacitance using atomic-layer deposition in 
a Josephson junction to prevent the junction from enter- 
ing phase diffusion. We compared the escape dynamics 
of the capacitively shunted and unshunted junctions. In 
the shunted junction, TA process is preserved at all mea- 
sured temperatures; for the unshunted junctions we ob- 
serve phase diffusion. The TA model yields a satisfactory 
discussion for our observations in the shunted junctions. 
It may be possible that, by engineering this shunt capac- 
itance for the experiments on the phase biased Cooper 
pair pumps [2|, we can not only improve the sensitiv- 
ity of the detector junction but also protect the circuit 
from high frequency noise. However, attention has to be 
paid on the fluctuations in the phase bias due to thermal 
effects. 
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